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Abstract: We define smoothed Wilson loop operators on a four dimensional lattice and 
check numerically that they have a finite and nontrivial continuum limit. The continuum 
operators maintain their character as unitary matrices and undergo a phase transition 
at infinite N reflected by the eigenvalue distribution closing a gap in its spectrum when 
the defining smooth loop is dilated from a small size to a large one. If this large N phase 
transition belongs to a solvable universality class one might be able to calculate analytically 
the string tension in terms of the perturbative A-parameter. This would be achieved by 
matching instanton results for small loops to the relevant large- iV-universal function which, 
in turn, would be matched for large loops to an effective string theory. Similarities between 
our findings and known analytical results in two dimensional space-time indicate that the 
phase transitions we found only affect the eigenvalue distribution, but the traces of finite 
powers of the Wilson loop operators stay smooth under scaling. 
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1. Introduction. 

Most of the concrete results of this paper were obtained by methods of lattice field theory. 
These results make up a first component of a larger, rather ambitious and quite speculative 
program which is defined in the continuum, independently of any regularization. The 
outline of the program makes up the bulk of this section. The intention is to put the 
presentation of the lattice results that follows in later sections in perspective. Throughout 
the paper, we shall get back to continuum language to emphasize what structures are, in 
our view, continuum properties. 
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Pure four dimensional SU(N) gauge theories are strongly interacting at large distances 
and weakly interacting at short scales. In SU(3) the cross-over between these two regimes 
is relatively narrow and there are no good analytic approximation methods that apply to 
it; nor do we have a continuum way to calculate at strong coupling. The cross-over narrows 
as the number of colors, N, increases. 

Here, in a direct continuation of the line of thought of some old work p], 0| 3 we report 
on an improved line of attack on this problem and take several steps in its direction. The 
old idea was that instantons indicate that in the limit of infinite N, at a particular sharply 
defined scale, there is a phase transition, and weak and strong coupling physics can be 
connected at that point. This idea has not been successfully implemented because too 
many parts were either wrong or missing. Much has been learned in the interim period 
and we would like to try again. Let us now briefly overview the new ingredients of our 
approach: 

While the guess that the cross-over regime collapses to a point and becomes a phase 
transition at infinite N is old, we understand now that the scale at which the transition 
occurs is not a universal quantity that applies to all observables simultaneously; nor does it 
necessarily affect any local observables: it is not seen in the free energy density for example. 
The most interesting large N transitions do not have to be "bulk" phase transitions. 

Recent lattice work on large N has produced such "bulk" transitions at global critical 
scales in the planar limit of QCD on a Euclidean torus, but the associated transitions take 
the system out of the zero temperature, infinite volume regime and depend on the global 
topology of space-time. This is not going to directly explain how a few particle high energy 
process hadronizes, for example. 

Recent lattice work has also indicated that specific non-local observables undergo large 
N phase transitions as they are dilated |3|, Q and this happens while the Euclidean space- 
time volume stays infinite in all directions. The critical scales are observable dependent; 
one can imagine associating with these observables an entire family of dilated images and 
somewhere along the dilation axis the large N transition takes place. 

If there were a string equivalent to planar QCD, one would guess that the dilation axis 
is somehow related to a fifth coordinate as has been often hypothesized |j, A meaningful 
higher dimensional formulation presupposes a local definition of the theory in the extended 
space; we shall have something to add to this speculation at the end of the paper, but it 
will be vague. 

Another new ingredient we add to the old program is the modern understanding of 
universal features of certain large N transitions in matrix models. Although, on the face of 
it, there is no concept of locality in the index space associated with the symmetry group, 
the models admit a 't Hooft double line expansion for two-index representations classified 
by the genus of two dimensional surfaces. A concept of two-dimensional locality applies to 
the physics on these surfaces and the associated Liouville type models represent universality 
classes of finite, but varying degrees of stability associated with the large iV transitions @. 
Most of these models are exactly solvable. 

We hope that the transitions we shall exhibit belong to some such universality class. 
We have a guess for the relevant universality class and it has an exactly solvable matrix 
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model representative. This means that at finite but large enough TV, when the transition 
dissolves into a cross-over, there is a description of well chosen observables undergoing the 
cross-over in terms of universal scaling functions with a dependence on a limited number of 
non-universal parameters that typically come in as boundary conditions on the equations 
determining the universal, known scaling functions. 

It is through these parameters that the cross-over regime connects to the weak and 
strong coupling regimes on either side. The four dimensional dilation parameter, in the 
vicinity of its critical value plays the role of a coupling that turns critical in the relevant 
version of matrix model description of the large N transition. By the mechanism of a (dou- 
ble) scaling limit the dilation parameter gets connected to a scale in the two dimensional 
physics that codifies the universal features of the large N transition. This connection will 
involve a nontrivial power and an associated dimensionful parameter. 

The weak and strong coupling regimes are governed by their own universal behavior, 
field theoretical on the weak side, and (we guess) effective string like on the strong coupling 
side. Thus, the free, nonuniversal, parameters needed by each universal description, field 
theoretic, large N and string-like are matched to each other as the operator is dilated, 
and several universality mechanisms coexist in the same system, controlling separate scale 
ranges. 

If the above scenario holds at N = oo, we can apply it to pure SU(3) gauge theory, and, 
hopefully provide a calculation of the string tension in terms of Aqcd which is accurate 
to order « 10 percent. For a physicist, this might be of more value even than a 
mathematical proof of confinement. We are far from realizing this goal. 

The first hurdle to overcome is to find good observables; we claim that the right 
question is not whether or not there is a large N phase transition in general, but rather, 
can one find good observables which undergo such a transition and can be controlled 
(by different theories) at very short and at very long distance scales respectively. Local 
observables are unlikely to undergo the transition. Therefore the cherished analyticity 
properties in four dimensional Euclidean momentum space stay protected, maintaining the 
basic unitarity and covariance properties of the field theoretical description of the physics. 

In this paper we define a class of smeared Wilson loop observables on the lattice with 
the intent that they have finite continuum limits and obey the following properties: 

• It is correct to think about the parallel transport round the loop as being realized 
by an SU(N) matrix. This statement has to hold after the renormalization of the 
theory and of the operator have been carried out. 

• When we scale the loop, a phase transition takes place in its eigenvalue distribution 
if N = oo. 

• The smeared loop operator is calculable in perturbation theory when it is scaled 
down. 

• Its shape dependence is well described by an effective string theory when it is scaled 
up. 
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• It is possible to see by continuum non-perturbative methods (e. g. instantons) that 
the phase transition is approached when the loop is dilated up from being very small, 
so one can hope to calculate the matching at the weak coupling end of the cross-over. 

At the large scale end one also needs a way to calculate a matching, and this is presumed to 
work by the large scale end of the cross-over having some string description (coming from 
the relevant universal matrix model) which can be joined smoothly to the right effective 
four dimensional effective string theory. That effective string theory could be, for example, 
similar to ||. This we have not yet established, and is left for future work. Nevertheless, 
on the basis of the results we do have, we believe to have found non-local observables which 
are "good" enough. 

In its essence the calculation we have outlined is an attempt to extend the logic of 
effective field theory |J calculations to the strong-weak cross-over in pure gauge QCD. 
Also, the program is in the same spirit as the old CDG approach to instantons fl(| where 
they tried to use instanton gas methods to set the parameters in an effective bag-model 
description of hadrons. 

If our scenario works for pure glue, it would be possible to re-introduce fermions at 
a later stage, at least as probes with no back reaction and see whether one can construct 
hadronic observables containing fermionic fields that also undergo large N transitions. 
In particular it would be interesting to find a similar calculational program for sponta- 
neous chiral symmetry breaking. One could check this numerically using exactly chiral 



fermions [11]. 



In the next section we introduce our observables on the lattice and present our nu- 
merical results in support of the claim that our observables have a continuum limit and 
undergo large N transitions under dilation. In subsequent sections we discuss a contin- 
uum view of the lattice construction and our guess for the large N universality class. The 
presentation becomes by necessity progressively more speculative. Ultimately, we hope to 
reduce the role of lattice simulations to a strictly qualitative one, providing support for the 
assumptions about the various universality classes that control the different regimes our 
calculation needs to connect. The calculation in itself would be entirely analytical. It goes 
without saying that any analytical step would need to be checked against numerical work, 
as this is unchartered territory. 



2. Overview of lattice study. 

We use an ordinary single plaquette Wilson action defined on a hypercubic lattice. We 



exploit twisted reduction [12] which allows a decent interval of lattice spacings to be realized 
on a small lattice of size V = 7 4 . For twisting to work best, we need N to be given by 
4 times a prime number. We work mostly at N = 44 with checks for finite N effects at 
N = 28. Our simulation method employs a combination of heat-bath and over-relaxation 
updates and equilibration is achieved in reasonable lengths of computer time. In most 
of our analysis we can ignore all statistical errors. As with any lattice work with few 
precedents to rely on, a substantial part of the effort is to identify useful ranges of the 
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many parameters that enter numerical work, so that the relevant physics questions can be 
addressed. Throughout, we use b for the lattice gauge coupling which is the inverse bare 
't Hooft coupling; relatively to the conventional lattice coupling (3, b already has the right 
power of N extracted: 

b = -A = -J— (2.1) 

2iV 2 g 2 YM N y ' 

For clarity purposes, we ignore a technicality in our presentation: For a twisted simulation 
b is taken as the negative of its value in the untwisted case; we shall quote the conventional, 
positive value of b. 

2.1 Smeared lattice Wilson loops. 

We consider only rectangular loops restricted to one plane of our hypercubic lattice and 
most of our work is done with square loops. However, a continuum approach would be 
concerned with the class of observables that contain smooth loops of arbitrary shape. 
Specifically, explicit calculations in the continuum would best be applied to circular loops 
contained in one flat plane in R , as they have maximal symmetry. 

In order to eliminate the perimeter and corner divergences [|l3| we construct the Wilson 
loop operators associated with our curves not from the SU(N) link matrices Unix) gener- 
ated by the probability distribution determined by the single plaquette Wilson action, but, 
rather, from smeared versions of those. We employ APE smearing |Q] , defined recursively 
from n = 0, where the smeared matrix is equal to U^{x). Let E^m)^ ^ denote the "staple" 

associated with the link (x; f) in terms of the entire set of ui n \y; /) matrices. One 
step in the recursion takes one from a set U^'{x\ f) to a set ujT +1 \x; /): 

X^ +1) (x;f) = (1 - \f\)Ujr\x; f) + 

U^Xx-f) = Af+ 1 )(*;/)^== = L=_= (2.2) 

V[X? +1) (x;/)]txi n+1) (x;/) 

In the simulation, one never encounters a situation where the unitary projection in the 
above equations stalls because is singular. In other words, smearing is well defined 
with probability one. For a twisted simulation / is taken as the negative of its value in the 
untwisted case; we shall quote the conventional, positive value of /. 

U^Xx; /) transforms under gauge transformations the same way as U^{x) does. Our 
smeared Wilson loop operators, W[Li, 1% /; n] are defined as ordered products round the 
L\ x L2 rectangle restricted to a plane. L a are integers and give the size of the loop in units 
of the lattice spacing. When the traversed link starts at site x = (x±, X2, x%, X4) x M G Z 
and connects to the neighboring site in the positive direction /i, x + n, the link matrix 
is U^ n \x; /), while when this oriented link is traversed in the opposite direction, the link 
matrix is U^ n \x; /). W depends on the place the loop was opened, but its eigenvalues do 
not. The set of eigenvalues is gauge invariant under the fundamental gauge transformation 
operating on U^x). 

We established by numerical means, to be described in greater detail later, that as 
L%, L2, /, n are varied, at specific lattice couplings, the spectrum of W[Li,L2, ; /; n] opens 
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a gap for very small and/or very smeared loops. This gap closes for very large and/or very 
lightly smeared loops. 



2.2 Lattice large N transitions survive in the continuum limit. 

Even restricting to just square loops, our lattice operators depend on many parameters: 
there is an overall scale, the side of the square, L; there is a gauge coupling b which fixes 
the lattice spacing in physical units, and therefore also controls the scale; there are two 
parameters, / and n that control the APE smearing and are responsible of the taming 
of the perimeter and corner divergences in the traces of the operator in all irreducible 
SU(N) representations simultaneously. Our construction implicitly adjusts the infinite 
number of undetermined finite parts associated with the renormalization of the trace in each 
representation so that the entire collection can be viewed as coming from one (fluctuating) 
SU(N) matrix. This ensures that backtracking segments of the loop make no contribution 
to the operator. 

Obviously, there is a substantial amount of arbitrariness in the choices we made when 
defining the regularized Wilson loop operators. If the construction is kept more general, 
this arbitrariness would reflect itself in an invariance under a large set of transformations 
in the continuum, perhaps under the form of a generalized renormalzation group equation. 

Our first objective is to find out how to adjust the parameter dependence in W such 
that N = oo transition points identified on the lattice, where the gap just opens, survive in 
the continuum limit in which the lattice coupling b is taken to infinity together with L\ i in 
such a way that the physical lengths l a = L a a(b) are kept fixed. a(b) is the lattice spacing 
at coupling b, and we express it in terms of the physical finite temperature deconfining 
transition t c , by 

a(b) = ^ (2.3) 

where is the lattice size of the compact direction on an S 1 x R 3 lattice corresponding 
to the finite temperature transition critical coupling b. 

We set the number of smearing steps n to be proportional to the perimeter square (we 
restricted the loop sizes to even Li + L 2 ), n = (LlH 4 L2)2 . The physical sizes of the loop are l a 
and L a = ^jjy. We have set n = (-M fe )+ L2 ( & )) because in physical terms the product fn is 
a length squared. This can be seen from the simple perturbative argument to be presented 
in the next subsection, or, more intuitively, by observing that smearing is a random walk 
that fattens the loop and the thickness grows as the square root of the number of smearing 
steps. Our choice for n makes / a dimensionless parameter in the physical sense; on the 
lattice / is actually bounded to an interval of order one. 

We analyzed four sequences of W[L±, L2; f;n = ( Ll + L2 ) ] keeping the physical size 
fixed at four values, one for each sequence. We set L\ = L2 = L. By various means, to be 
described later, we varied / and found the critical value of / at which each member in each 
sequence just opens a gap in its spectrum at eigenvalue —1. (So long as CP invariance is 
maintained, for any loop, the eigenvalue probability is extremal at ±1.) The sets of critical 
values for each sequence were then extrapolated to the b — ► 00 limit. In each sequence 
the physical size of the loop is kept fixed by definition. This gave four continuum critical 
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values of / for the four different physical loop sizes we selected. Throughout, the values of 
/ used were within reasonable ranges, well within the maximally allowed lattice range of 
< / < 0.75 [15|. That it is actually possible to arrange that only such reasonable values 
of / are needed in the vicinity of the large TV transitions is our most basic result and a 
nontrivial feature of the dynamics of planar pure Yang Mills theory. The collection of the 
continuum limits of the critical values of / at the transition points associated with each 
sequence identified by its physical scale defines a critical line, / c (0- For any finite b there 
were finite lattice spacing effects, but they were sufficiently small to keep all the lattice 
approximations to / in a good range. 

Thus, we conclude that we have found a continuum critical line f c {l) in planar QCD 
where our continuum Wilson loop operators have a well defined spectrum and that spec- 
trum is critical in the sense that its eigenvalue density, p{9) is nonzero for any 6 ^ ir but 
p(ir) = 0. Here, we denote the eigenvalues of the Wilson loop operator as e tS , where 9 is 
an angle on a circle. 

It is not necessary to keep n quadratic in L a for Al a >> 1; the loops that have a 
critical / are of sizes I for which Al is of order unity. For much larger loops, fn should be 
allowed to become proportional to a physical scale of order and not increase further as 
the physical loop size goes to infinity. All that is important is that for loop sizes of order 
j, n be adjusted as above, and that for small loops we let fn become small too. For small 
loops we wish to enter a perturbative regime, and therefore we should not eliminate the 
ultraviolet fluctuations at short distances, where they dominate. 

Once the existence of the f c {l) line has been established in the continuum limit, with 
^ > 0, we can be reasonably sure that dilating a W[L, L; /; n] so that at lattice coupling b 
it passes through a point L = / = f c (l) + 0(a 2 (b)),n = L 2 the spectrum of this W will 
go through the transition. We still need to make sure that once the gap is closed, the density 
at —1 starts building up continuously from zero; this would mean that the transition is 
continuous and likely to possess a universal character. If the transition were discontinuous, 
like, for example, the finite temperature transition in Polyakov loop operators, it is unlikely 
that a scaling regime exists on either side of the transition and a single scaling regime 
connecting the two sides (of the type we are hoping for) is ruled out. This is obviously the 
case in the finite temperature case because the transition there is not induced by taking 
TV" — > oo, but exists for all N, and starts being first order at N = 3. 

We addressed this point by carrying out additional simulations along another line in 
the physical (/, I) plane, which is not a constant I but, rather, has / varying by a substantial 
amount. We picked a convenient interval of lengths I and a line in the plane (/, I) given by 
/ = — fil 2 with constant f\ > chosen so that / varies between f\ and fi with f\ 0.1 
and /2 ~ 0.2. This line intersects / c (/) at an angle not too close to zero, thus enhancing 
the variability in the spectrum as one takes W along the line between the two extremal 
lengths. There is no particular physical meaning attached to this second line, except that 
along it the loop is dilated and in addition, the amount of smearing decreased; both effects 
help closing the spectral gap that is present when the loop is small. 

To go along the continuum line we fixed the lattice loop-size L and varied b. We kept 
n = L 2 also fixed, but varied / according to the formula / = f(b) = fo — f\(La(b)) 2 . 
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Following W along this line we established that indeed the spectrum of W opens a gap as 
the band of lattice curves converging to f c (l) is crossed. 

We avoided working close to the transition and estimated independently on each side 
where the corresponding phase would end, using various extrapolation methods: On the 
side where there is a gap, we studied the decrease of the gap to estimate when it will close. 
On the gap-less side, we estimated when p(n) would vanish. 

We repeated the procedure for three increasing values of L. In each case one covers 
similar (overlapping, but not identical) ranges of physical size I. In the common portions 
of the range, the larger the lattice value of L is, the smaller the lattice effects ought to 
be. We found general consistency with a scaling violation that goes as the lattice spacing 
squared. 

The separate extrapolations in each phase produced independent estimates for end 
points of each phase which came out reasonably close to each other. This is consistent 
with a continuous large N transition at which the eigenvalue distribution of W opens a 
gap. The determinations of the critical points from each side were ordered in a way that 
was inconsistent with a discontinuous transition in which p(ir) decreases to a non-zero value 
and then drops to zero, jumping into a phase with a finite gap around tt. Thus, we conclude 
that the transition is continuous and one can hope for some universal description in the 
vicinity of the transition point. Actually, it is hard to imagine what kind of dynamics could 
have produced a first order large N transition in the spectrum. 



2.3 Smearing removes ultra-violet divergences. 



The effect of smearing is easy to understand in perturbation theory where one supposes 
that each individual step in the smearing iteration can be linearized. Writing uf?\x\ /) = 
ex.p(iA^ \x; /)), and expanding in one finds [15], in lattice Fourier space: 



(«;/) = EV(?)4 n) fe/) (2.4) 

V 

with 

M9) = /(?)(V-% + ¥ ( 2 - 5 ) 



where = 2sin(^ i ) and 



f(q) = i - 4~ 2 (2.6) 

o 



The iteration is solved by replacing f(q) by f n (q), where, for small enough /, 

r(q) ~ e"^ 2 (2.7) 



The condition < f(q) < 1 for all momenta, requires < / < 0.75. To get some feel 
for the numbers, a 4 x 4 loop, at a value of b that makes it critical, will have n = 16 and 
/ ~ 0.15, which makes ~ 0.4, which is an amount of fattening over a distance as large 
as the loop itself. This is why much larger loops should not be smeared with an fn factor 
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that keeps on growing as the perimeter squared; rather, for a square loop of side L, for 
example, the following choice would be appropriate: 



Here, M is a physical mass m expressed in lattice units, M = ma(b). m is a typical 
hadronic scale chosen so that at the large N transition we found, ML is less then 0.01, say. 

To reassure ourselves that such an extension of the cutoff dependence of / indeed 
would connect smoothly to large loops which show confinement with the known value of 
the string tension, we roughly estimated the string tension from our smeared loops, at the 
large size end of the ranges we studied. We found that the string tension comes out close to 
that of bare lattice Wilson loops, and this shows that our smeared Wilson loop operators 
are acceptable pure SU (N) gauge theory string candidates. 

In summary, the effect of smearing on the propagator of Ap(x; f) is to eliminate the 
ultraviolet singularity that gives rise to perimeter and corner divergences in traces of bare 
Wilson loop operators. This works when the product fn varies as a physical length squared 
when the continuum limit is approached. For small loops we choose this physical length 
to be proportional to the loop itself. For large loops the physical length stabilizes at some 
hadronic scale. In this way, our operators become sensitive to the short distance dynamics 
when the scale is very small, but, when confinement becomes the dominating effect, for 
very large loops, the string tension of the smeared Wilson loops is the true, universal, 
string tension of the theory, unaffected by the smearing. The leading term in the effective 
string theory, describing the shape dependence of the Wilson loops when the loop is very 
large, is expected to be universal and the string tension enters only in a trivial manner, as 
a dimensional unit. 

We stress that we are not executing a renormalization group transformation; our Wil- 
son loop operators are defined for infinitely thin loops in space-time and for all loop sizes; 
rather, we are applying a regularization whose objective is to make the Wilson loop opera- 
tors finite in the continuum. If we dealt with each irreducible representation separately, the 
elimination of perimeter and corner divergences would leave behind an infinite number of 
arbitrary finite parts. This freedom has been implicitly exploited by imposing the require- 
ment of unitarity of the Wilson loop operators, which ensures what Polyakov has termed 
the "zigzag" symmetry ||. Preserving the unitary character of the Wilson loop operator 
in turn makes the concept of an eigenvalue distribution well defined in the continuum and 
sets the stage for possible large N phase transitions. The zigzag symmetry also influences 
the form of the lattice loop equations at infinite N, which are obeyed by traces of bare 
Wilson loops in the fundamental representation; the effective string theory for large loops 
needs not obey the zigzag symmetry, but this symmetry could provide a useful constraint 
on its couplings. 

3. Numerical details of lattice study. 

In the previous section we described what our strategy was to address the question of 
existence of a continuous limit for the operators and the large N transition point. Also, we 



/ 
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made an effort to convince ourselves that the large N transitions indeed are continuous. 
We now proceed to describe the exact details of our simulations and the numbers that 
came out. 

3.1 The determination of f c (l). 

3.1.1 Estimating the gap at infinite N. 

At infinite N, small W loops will have a spectral gap on the lattice. We need a method 
to extract that gap from numerical data at finite N. As we dilate the loop, and/or reduce 
the amount of smearing, the gap shrinks. We need a way to extrapolate to the point where 
the gap vanishes. We wish to keep away from measurements for too small gaps, and enter 
that regime only by extrapolation. 

Getting the infinite N value of the gap is difficult because subleading effects are large 
at the gap's edge. It is reasonable to assume that the universal behavior at the gap's 
edge is the same as at the edge of a gap in simple hermitian matrix models. Under this 
assumption, we have separate predictions for the behavior of the eigenvalue density and 
for the extremal eigenvalue. One can check numerically how well the universal forms work 
and we find them to work quite well. To fit onto the universal behavior we need to adjust a 
parameter that corresponds to the gap size at infinite N. This provides the determination 
of the gap sizes we use subsequently. 

It is only on the side of the transition where the eigenvalue density has a gap that there 
is a regime where large N corrections are dominated by universal functions; therefore, we 
are able to get good estimates for the eigenvalue gap size at infinite TV from data obtained 
at one single, but large N. For good measure, we checked that the predicted rational power 
dependence on N holds in the regimes in which we did our analysis. A comparison with 
some simulations at N = 28 convinced us that the power dependence on N is already at 
its known universal value at N = 28, 44. 

We studied the eigenvalue distribution of W[L, L; f;n = L 2 ;N] on the lattice for 
several values of L and b at a fixed large N. Let e tdk , k = 1, • • • , N be the N eigenvalues 
of a fixed L x L loop on the lattice, — ir < Q\ < 62 < ■ ■ ■ < 0jv-i < @N < tt- The sum of all 
angles is zero mod 2n. On a lattice of V sites there are 6V different L x L Wilson loops 
with a fixed orientation; we computed the eigenvalues of all such loops. In this manner, 
we obtain a histogram of the distribution pf(9;U) of all eigenvalues in a fixed gauge field 
background U. We also collected data for a histogram of the distribution p m (0N;U) for 
the largest eigenvalue in a fixed gauge field background. In addition, we record the mean, 
variance, skewness and kurtosis pertaining to the largest eigenvalue: 




K N (U) 



S N (U) 



(e%) u -3(e 2 N ) u (e N ) u + 2(6 N )f J 

(e%)u - ±{e%)u{e N )u + Q{e 2 N )u{e N )l - 3(0^ 



- 3. 



(3.1) 
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If the eigenvalue distribution of W has a sizable gap, it is reasonable to ignore the 
influence of the two edges on each other. Then, the universal behavior at the gap's edge 
ought to be the same as that of the edge of the gap in simple Gaussian hermitian matrix 
model. Under this assumption, we compared our data with the behavior of the extremal 
eigenvalue 16] and the eigenvalue density close to the edge JT7|. 



The universal distribution of the extremal eigenvalue, p(s), is given by 



pO) 



OO 



2 



x)dx 



e - jr( x ~ s )i 2 ( x ) dx (3.2) 



with q being the solution to the Painleve II equation, 



fq 
ds 



- 1 ^=sq + 2q 3 ; (3.3) 



satisfying the asymptotic condition 

q(s) ~ Ai(s) as s — > oo. (3.4) 

The mean, variance, skewness and kurtosis of the above distribution are —1.7710868074, 
0.8131947928, 0.2240842036 and 0.0934480876. In order to match the lattice data for the 
distribution of On to the above universal form, we define s according to 

On = - + tt(1 - g) (3.5) 
a 

and set 



/0.8131947928 



a 



V a N 

l r 1.7710868074. , . 

g = l--\ptN + ] (3.6) 

7r a 

a is a non-universal scale factor and g is half the arc length of the gap on in the eigenvalue 
distribution on the unit circle, centered at —1. The quantities, a N and hn are the averages 
over the gauge field ensemble of the previously defined quantities, a N {U) and hn(U). 
From the lattice data for a N (U) and /j,n{U), we obtain an estimate for a and g with errors 
extracted by the jackknife method. Using these estimates, we compared the scaled lattice 
data with p(s) for fixed /, b, N and V. One such comparison, at / = 0.15, b = 0.361653, 
L = 4, N = 44 and V = 7 4 with a gap of 0.1229(4), is shown in Figure |] and there is 
good agreement. In spite of the fact that data is needed over the entire range of s to get 
a good estimate of the skewness and kurtosis, we find agreement between lattice data and 
the universal estimates to within 5% for the skewness and to within 20% for the kurtosis 
for a wide range of gap values. The estimate of the skewness and kurtosis do not depend 
on the extracted values of a and g. 

The main result of the above numerical analysis is the half gap size g, for fixed /, L, b, 

g(f,L,b). 
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Figure 1: A comparison of the distribution of scaled extremal eigenvalue obtained on the lattice 
(circles) at / = 0.15, b — 0.361653, L — 4, N — 44 and V — 7 4 with that of the universal distribution 
(solid curve). 



There is another way to obtain an estimate for g: A recent result provides the asymp- 
totic corrections to the eigenvalue density of Gaussian hermitian matrices [|l7|] . Using that 
formula, we fitted he eigenvalue density in terms of the scaled variable in ( |3.5| ) to 



p f (s) = c { [Ai'(s)} 2 - s[Ai(s)} 2 } + d{3s 2 [Ai(s)} 2 - 2s{Ai'(s)} 2 - 3Ai(s)Ai' (s)} (3.7) 



The variable c is fixed by the overall normalization and 



l 

" 2(W 2 /3 



in the asymptotic 

formula [17] where iV is the size of the matrix. The parameter d could be non-universal; also 
c could in some sense be non-universal, since it depends on the total number of eigenvalues, 
and it is not clear whether the Gaussian hermitian matrix model matches onto our model 
with the same N in both cases. Therefore, we fitted the lattice distribution to the above 
asymptotic formula with c and d as independent variables in addition to the two variables 
appearing in ( |3.5| ), producing a second determination of g(f, L, b). 

Figure [2] shows the comparison of the full eigenvalues distribution with the universal 
distribution near the edge. The universal distribution captures the tail for the same reason 
that the distribution of the extremal eigenvalue is in agreement with p(s). But we also see 
good agreement with the first two oscillations near the edge and this includes more than 
the extremal eigenvalue. In general we found the determination of the gap to be consistent 
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Figure 2: A comparison of the distribution of eigenvalue distribution obtained on the lattice 
(circles with dashed lines) at / = 0.15, b = 0.361653, L = 4, N = 44 and V = 7 4 with that of the 
universal distribution near the edge (solid curve). 

with the one obtained using only the extremal eigenvalue. For an error estimate, we only 
used the extremal eigenvalue method. 

3.1.2 Extrapolating to zero gap. 

The next step is to extract the critical value of / at a fixed b and L at infinite N. This is 
done by assuming a fit of the form 

g(f, L, b) = A(L, bWf-F c (L,b) (3.8) 

The square root dependence on external parameters that enter the probability law in an 
analytic manner is generic in matrix models and well supported by the data. 

A plot of g 2 (f,L,b) as function a / at b = 0.361653, N = 44 is shown in figure ||[ 
F c (L,b) is extracted using a range of / such that the gap, g(f,L,b), is not too close to 
nor is it too far from 0. If we go to too small gaps we fear unaccounted for rounding effects; 
if the gap is too large, there is little reason to use the universal square root behavior we 
postulated. 

Figure |3] also shows the fit along with the estimate for F C (L, b). 
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Figure 3: Plot of g 2 (f, L, b) as a function of / at b = 0.361653, N = 44 and V = 7 4 . A fit to the 
data along with the extrapolated value for F C (L, b) is also shown. 

3.1.3 The continuum limit of F c (L,b). 

We now proceed to take the continuum limit of F C (L, b). We used two loop tadpole improved 
perturbation theory to keep the physical loop fixed as we varied b. Using we have a 
numerical formula for the critical deconfinement temperature in lattice units, in terms of 
tadpole improved perturbation theory, with e(b) being the average value of the plaquette 
for untwisted single plaquette gauge action at coupling b. 



' 0.26 



T c (b) 



11 



48vr 2 6/ 



51 



121 24^ b, 

e - n— ; b! = be{b) (3.9) 



This will be used to determine the lattice spacing in physical units. 

We set lt c = LT c (b) and accumulate data at several values of the physical loop size: 
lt c = 0.42,0.48,0.58,0.64. At each of these physical sizes, we went through a sequence of 
three to four lattice spacings a(b), and, for each b, we obtained a value for F C (L = b), 
using the methods explained above. 

For each lt c , the F C (L = ^4y,&) values were then extrapolated to the continuum as- 
suming 

F C (L = -L,b) = f c (lt c ) + 0(T C 2 (6)). (3.10) 
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Figure 4: Plot of F C (L — b) as a function of T c (b) for four different values of lt c . Also shown 
are fits to the data producing continuum values for f c (lt c ). 



where we recall that T c (b) = a(b)t c . (We shall use the slightly inconsistent notations 
f c (J) and f c {lt c ) for the same quantity, where the argument is either dimensional or made 
dimensionless by the appropriate power of the physical deconfining temperature t c .) 
The extrapolated results for f c (lt c ) are shown in figure ||. 

Our final results are plotted in the (/, I) plane. The critical line f c (l) = fc(ltc) divides 
this plane into two as shown in figure [|. In the upper part Wilson loops have a gap in their 
eigenvalue distribution, whereas in the lower part there is no gap. The band bounded from 
below by the critical zero-lattice-spacing line indicates the range of finite-lattice-spacing 
effects. 



3.2 A study of both sides of the transition under scaling. 

Finite lattice size effects are evident in figure |5[ It is impractical to just scale the loop 
because the curve f c {l) varies when I is changed in the appropriate range by less than the 
lattice corrections. To see the effect of scaling we need to go on a line in the (/, I) plane 
that is not horizontal, which means we are going to vary the smearing as a function of the 
physical loop size. As usual, we trade I for the dimensionless variable lt c . 
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Figure 5: The critical line in the (f,lt c ) plane. Also shown is the crossing line used to investigate 
the phase transition induced by scaling the Wilson loops. 

We chose the line (we shall refer to this line as the "crossing line") 

/ = 0.25 - (0.6128/t c ) 2 (3.11) 

shown in figure ^ as it cuts across the entire finite band of lattice corrections bounded by 
the critical line f c (lt c )- 

Consider a fixed loop of size L on the lattice for various lattice couplings b such that 
f(b) is given by fljUlj ) with lt c = LT c (b) and T c (b) given by ( $S§ . We wish to study the 



nature of the transition in the eigenvalue distribution along this line. 

We already know that the distribution should have no gap for I > l c whereas it should 
have a gap for I < l c . In addition to getting an estimate for l c t c , we can also investigate 
the universality class of this transition by keeping L fixed and going along the crossing line 
given by (3.11). We set L = 3,4,5 and worked in the following ranges of b: for L=3, in 



0.345 < b < 0.365, for L=4, in 0.346 < b < 0.368 and for L=5, in 0.351 < b < 0.369. In 
each range, as b decreases the physical loop is scaled up. In this way one gets overlaps at 
the same physical loop size, once represented by a smaller value of L and a smaller value of 
b, and other times, by correspondingly larger values of these parameters. Because of scaling 
violations one will not get identical eigenvalue distributions; we need to see that, as b and 
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L increase, these corrections decrease and that a finite limiting eigenvalue distribution is 
approached. This is the eigenvalue distribution associated with the continuum loop of size 
I. 

In Figure [5] we display, for each value of L = 3, 4, 5, the two critical lengths one obtains 
when extrapolating from each phase independently as triangles on the crossing line. As L 
is increased, we see that the distance between the two members of each pair decreases as 
L increases. Also, the middle points of the segments connecting the points in each pair 
appear to converge as L increases. This indicates that in the continuum limit there will be 
one single transition point on this line, a loop size at which a continuous transition in the 
large N spectrum of W occurs. 

3.2.1 Estimating the critical size from the small loop side. 

We again obtained estimates for the infinite N gap, using the extremal eigenvalue distri- 
bution (p(s)) and the eigenvalue density (j>f(s)) near the edge. The estimates for the gap 
came out consistent with each other and the results are plotted in figure |6[ The critical size 
of the loop is extracted by again assuming a square root singularity written in the form 

g 2 (lt c )=Alnj- (3.12) 

For each value of L = 3, 4, 5 we got an estimate for l c t c which is shown in figure [| One 
sees that the lines L = 4, 5 are closer to each other than the lines L = 3, 4; the ratio of 
the approximate distances between the two pairs of lines is in rough agreement with a 
hypothesis of order a 2 (b) corrections. 

3.2.2 Estimating the critical size from the large loop side. 

For operators with no gap we need a method to extract the infinite N value of the eigenvalue 
density at angle it, p(ir). Shrinking the loop and/or smearing it more makes p(ir) decrease. 
We also need a method to extrapolate to the point where p(ir) vanishes, without getting 
too close to it. 

The estimation of the infinite N value of p(jr) also suffers from large subleading effects, 
of order -i typically. There is a rapid fluctuation in the eigenvalue density reflective of the 
repulsion between eigenvalues. The rapid fluctuation is eliminated by Fourier transforming 
the numerically obtained binned eigenvalue distribution with respect to the angle for wave- 
lengths kept large relative to -k • This is equivalent to evaluating the traces of W k for 
consecutive values of k, all kept much smaller than N. Keeping only those Fourier modes, 
we get our estimates for p(n); checks against a simulation at lower TV's shows that we are 
determining p(ir) to reasonable accuracy. 

The extrapolation of p(ir) to zero as the loop size is decreased and/or its amount of 
smearing is increased is done in a similar way as when we deal with the gap. We find that 
the leading behavior is again square-root like, but this holds in a small range which can 
be extended by adding more powers of (I — l c ). The determination of l c by this method is 
seen to be quite robust. 
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Figure 6: The square of the gap as a function of ln(7i c ) along the crossing line. 



In two dimensions the large N eigenvalue distribution for loops of arbitrary size is 



known in the continuum [18]. As we shall elaborate later, it is plausible to assume that 
the behavior close to the large N transition is universal, and holds also for our p{6) in four 
dimensions. Therefore, we postulate that p(iv) vanishes as a square root at the critical value 
of any parameter that enters the underlying distribution controlling the smeared Wilson 
loop operators in an analytic manner. This is why we took p(n) to vanish as \Jl — l c as 
one approaches the appropriate critical size, l c , from above. 

As indicated above, for I large enough with respect to l c we extract the nonzero value 
of p(ir) by first representing the distribution for all angles as: 

N 

/5(0) = j]/jfecos(fc0) (3.13) 

fc=0 

Next, we truncate the Fourier series and use the truncated series to provide the value of 
p(7r). 

km 

p(ir) = (3.14) 

k=0 

The value of k m is chosen to be slightly less than N/2; this eliminates oscillations of 
frequency N. We checked that our choice for k m was reasonable by looking for stability as 
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Figure 7: The fourier coefficients fk as a function of fc for two different loops in the gapless region. 
The measurements were obtained using an L = 5 loop at N = 44 at two different lattice spacings. 

we vary k m . 

We found that for low enough &'s fk oscillates in sign from k = 1 onwards. As Figure 
shows, for k up to about ^- the coefficients are smooth in k, and the sign oscillation is 
obeyed. We shall come back to this property later. 

The values of p(w) we obtain are fitted to a square root in the critical region and the 
results are shown in figure |8| The estimates for l c t c from p(jr) tend to be higher than the 
ones coming from fitting the gap. This difference decreases as the spacing goes to zero. 

The critical values for l c t c obtained from the gap and p(ir) are shown in figure H One 
can see that the difference between the two estimates shrinks and one can also see that the 
average of the two estimates approach a finite limit as the lattice spacing decreases. 

The size and behavior of the scaling violations we see are consistent with usual expec- 
tations, reflecting contributions of higher dimensional continuum operators to the action 
and the operators. We leave for the future a more detailed investigation of these scaling 
violations, but note that the orders of magnitude are typical. 

3.3 General features of f c {l). 

In summary, using all of our data, a rough estimate for the critical size on our crossing 
line in the (/, I) plane, is l c t c 0.6. l c « y/a, where a is the string tension and we used 



- 19 - 



0.2 



0.1 - 



• N=44, 3X3 

1.33 ln(lt c /0.560) 

- N=28, 3X3 

-- 1.191n(lt c /0.579) 

♦ N=44, 4X4 

._. 1.601n(lt c /0.601) 

a N=44, 5X5 

.... 1.71 ln(lt c /0.626) 



i 1 1 1 1 1 1 — 71 — 71 — r~i r — I 1 1 1 r 




ln(lt c ) 



Figure 8: The square of the p(ir) as a function of ln(Z£ c ) along the crossing line. 



the approximate relation 0.6a = t c . In QCD units, l c is about half a fermi, in general 
accordance with the accumulating body of numerical data which sets the transition from 
perturbative to string-like behavior at rather short scales, of less than one fermi. Clearly, 
varying the crossing line, the loop shape and/or the details of the continuum construction 
of the Wilson loop operator will produce slightly different values of l c . However, the small 
variation over scales in the critical line f c (l) we find indicates that the transitions will 
occur in a restricted range of scales. Numerically, it seems that f c {l) asymptotes to a 
positive constant as / — ► 0. This excludes choices of Wilson loop definitions with l c t c « 1. 
However, if it is somehow determined by further simulations on very fine lattices that 
/c(0 ~~ * as I — > 0, such choices must be recognized as possible; in that case it is unlikely 
that an effective string description will hold for I too close to l c . Even if it were possible, 
we would not choose such a definition of regulated continuum Wilson loop operators. 

We have seen the scaling violations in the F c (L,b) approximations to f c (l) in a range 
of physical scales close to the large N transition point. It is impossible to decrease the 
physical size of the loops much more, because that would require too large values of b for 
the large N twisted reduction trick to work on a 7 4 lattice. Similarly, it is impossible to 
get to much larger loops because working at too small values of b takes the lattice system 
through a first order phase transition to a strong coupling phase which has no relevance to 
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continuum physics. 

We found that f c {l) flattened out for small loops. It is possible that /(/) does decrease 
to zero as I goes to zero, but if this at all happens, one would expect the effect to be 
logarithmic (i.e. f(l) ~ i log ^ i), since /(/) ought to behave roughly as the gauge coupling 
constant at scale /, which enters the perimeter divergence smearing is intended to eliminate. 
When I increases, / c (/) starts increasing more rapidly and eventually, f c {l) would violate its 
lattice bound. It is a strain to arrange for too large loops to have a spectral gap, because 
for large loops confinement is reflected by an almost uniform eigenvalue density. 

Note however that confinement is not necessary in order to have our large N transitions; 
even if the eigenvalue density does not approach a uniform distribution for very large loops 
it remains possible that large loops have no gap in their spectrum. Thus, a similar analysis 
might apply also to models that do not confine. 



4. The continuum view. 



We would like to formulate everything we have been doing in continuum language because 
the basic phenomenon is, we claim, not a lattice artifact. So, we need to understand in 
continuum language what smearing means, what the smeared observables are, and how one 
could calculate them using perturbative and non-perturbative continuum methods. Next, 
we would need at least to consider how the large ./V transitions would fit into a continuum 
description and whether they are at all acceptable in a unitary field theory. 

4.1 Smeared continuum Wilson loops and hints of a fifth dimension. 

Our first task is to understand what smearing means in the continuum. The answer is that 
the quantity fn plays the role of the fifth "time", t, in the Langevin equation when the 
noise has been set to zero: 

This can be seen by comparing equations (2.4-2.7) to (4.1) for small fn and small r. 



If we add a Zwanziger term |1£] to the right hand side (it has no effect on our Wilson 
loops or on any other gauge invariant observable) and move it to the left hand side, the 
equation of motion gets a five dimensional form, with x§ = r: 



F 5 „ = D^ u (4.2) 

When a noise term is added to the right hand side one gets the Parisi-Wu stochastic 
quantization of gauge theories [pC| ]. This setup has been shown to correspond to a topolog- 
ical, gauge invariant five dimensional theory theory |21], in which the fifth direction plays a 
special role (for example its units are length square as opposed to the units of the other four 
coordinates). There is Euclidean rotational invariance only in four dimensions. The bulk 
theory is trivial in the sense of a topological field theory, but the boundary has nontrivial 
dynamics. This setup is reminiscent of the concept of holography. The noise term is an 
essential ingredient in the construction, and it is an open and perhaps interesting question 
whether our work can be made to have a more precise relation to the one of pi] ]. 
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In our paper we only looked at loops at one common fixed r. This r determined the 
thickness of the loop which we took as uniform. But, it makes perfect sense to also consider 
loop operators where r varies round the loop, describing a string of varying thickness. 
Moreover, one would have to consider such loops if one looks for a loop equation that is 
local all around the loop. Indeed, the ordinary four dimensional equation of motion enters 
the r evolution and therefore it is possible that at the expense of considering loops of every 
possible thickness simultaneously one can restore the locality of the Migdal-Makeenko 
equation that we lost when we simply fattened the loops staying in four dimensions. 

In turn the Migdal-Makeenko equations are perhaps the right starting point to get at 
the fundamental (infinitely thin) large N QCD string, which would live now in five warped 
dimensions and might provide a description of four dimensional physics that is equivalent 
to that provided by four dimensional fat strings. This speculation suggests that our simple 
objective to get well defined continuum Wilson loops as unitary matrices in four dimensions 
did not lead us towards and extra dimension by sheer coincidence. 

4.2 Instantons close the gap. 

For a small loop the parallel transporter round it ought to be close to the identity matrix 
and have no eigenvalues around —1. At infinite N, one expects a total suppression of 
eigenvalue far from 1, and a large gap. In matrix models, for finite but large N, this 
suppression is typically exponentially small in N. 

Consider a circle of radius R in a plane and an instanton with center placed somewhere 
on the circumference of the circle. We can easily calculate the eigenvalues of the Wilson 
loop operator for the circle, opening it at the instanton center. 

The instanton gauge field in a regular gauge is given by: 

where 

U= X * + *-*. (4.4) 



Defining our circle as (x± — R) + x% = R , with x\p, = 0, we find: 

p e 'MA = e Mr(1 ~7?fi5ff) ' 3 (4.5) 

So long as the instanton is small enough and close enough to the periphery of the circle 
it would create some eigenvalue density at —1. Large instantons are irrelevant for AR « 1, 
so one does not need to face the well known infrared problem of instanton calculus if the 
loop is small. This is a e ~ N ' Const effect, likely to join smoothly onto the asymptotic behavior 
we would get from any matrix-model universal (double-) scaling function. This might be 
enough to select a unique solution to the differential equation of the matrix model and 
the latter may take us all the way to large R, where an effective string theory has an 
overlapping regime with the matrix model. Clearly, without some explicit calculations this 
is just a suggestion. 

Note that the continuum version of smearing, as we defined it, will not affect the 
instanton since it sets to zero the right hand side of ( f4.1| ). 
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4.3 A guess for the large N universality class. 

There are at least two obvious candidates for the universality class of the large N transitions 
in our smeared Wilson loops. If one just asks for the generic model dealing with a unitary 
matrix that opens a gap in its spectrum and which obeys a reality condition reflecting 
CP invariance, one will end up with a single unitary matrix model and the famous Gross- 
Witten |g] large N transition. The other immediate option is two dimensional QCD defined 
on a sphere of dimensionless area A which has the Douglas-Kazakov transition. The latter 
transition also occurs for two dimensional QCD on the infinite plane, but not for the 
partition function this time, but, rather for Wilson loop operators of area A; the transition 
on the sphere in one case and in observables in the other are mathematically identical. We 
now argue that the second option is preferred, both theoretically and numerically. 

It may be a bit confusing that we distinguish these two cases, as both are relevant 
to two dimensional QCD. The Gross- Witten transition is relevant to lattice QCD with a 
single plaquette action of the Wilson type. The single unitary matrix in the model is the 
parallel transporter round an elementary plaquette. As a function of the lattice coupling, 
this matrix opens a gap at some critical value as b is increased from zero and the trace of 
this matrix will be nonanalytic at the transition point. Now, consider the lattice model for 
which the analysis of Gross and Witten applies, but consider a loop larger than lxl. It 
is still true that its trace will be nonanalytic at the Gross- Witten transition point. But, 
the parallel transporter round the larger loop will open a spectral gap at a large value of 
6, and at that value the trace of the larger loop does not exhibit a nonanalytic behavior. 
In short, the nonanalytic behavior of the trace of the plaquette is a lattice effect which 
has no relevance to the continuum limit. Once we agree that we should avoid looking at 
lattice effects, it is obvious we should rather consider continuum two dimensional QCD as 
providing a representative of the large N universality class relevant to the transition we 
found for our smeared four dimensional Wilson loops W. 

In continuum two dimensional planar QCD the spectrum of Wilson loop operators, 
for smooth simple loops, does not depend on the the shape of the curve but only on the 
enclosed area. As the loop is scaled, all that matters is that the area changes, and at 
infinite N it is known that its eigenvalue distribution undergoes a transition at a specific 
area (measured in units of the two dimensional gauge coupling constant) where it just 
opens a gap at -1. However, the traces of any finite power of the Wilson loop operator stay 
analytic as the area goes through the critical point of the spectrum. The expectation value 
of an n- wound Wilson loop in two dimensions still depends just on its fundamental area, 
which we make dimensionless using the 't Hooft coupling, and denote by A. The explicit 
formula at large N is known [23|, given by the product of a polynomial and an exponential: 

±(TrW n ) = U^U2An)e- An (4.6) 

The dependence on scale hardly could be more analytic. 

The Gross- Witten universality class and the one associated with continuum QCD are 
related but not necessarily identical. Still, both candidates produce a Painleve II equation 
in the double scaling limit [24, |25|| , 
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Another point to keep in mind is that the string theory associated with the universal 
class of the Gross- Witten transition is quite degenerate and appears to differ significantly 
from any putative string theory dual of QCD [26]. On the other hand, there does exist 
a more tangible string description of two dimensional QCD due to Gross and Taylor [27] 
which also has sigma model representations [28|. The large N transition in Wilson loops of 
continuum two dimensional QCD on infinite space-time is of the same type as the Douglas- 
Kazakov transition in the partition function of planar two dimensional QCD defined on a 
sphere of area A [22]. 

There is another feature in favor of two dimensional QCD as the universality class 
and it has to do with instantons. When one takes the nonabelian instanton solution that 
we showed makes a non-zero contribution to the eigenvalue density at —1 for small loops 
in four dimensions and looks only at a two dimensional slice one gets two dimensional 
instantons that were argued to be related to the Douglas-Kazakov transition in [25|. 

All this adds up to putting the Gross- Witten transition in disfavor as a representative 
of the universality class of our transitions, and leaves us with the preferred option of two 
dimensional planar QCD. An explicit formula for the eigenvalue distribution of Wilson 
loops can be found in |18| and it is possible |2^] to connect this formula to ( |4.6| ). 

To see that we are not completely off in our guess we show a plot of several eigenvalue 
distributions we measure versus the exact infinite N continuum solution of Durhuus and 
Olesen @: 

Figure || shows shows how the critical size is traversed and one sees directly the effect 
of the transition and that the Durhuus-Olesen distribution seems to work pretty well on 
both sides of the transition. An attempt to fit to the Gross- Witten distribution would fail, 
since on the gap-less side the Gross- Witten distribution would always cross the value \ in 
Figure | at 9 = ±§. 

Another point in favor of the universality class is evident from equation 4.6: The 
prefactor has oscillating signs for large areas. These are nothing but the oscillating signs 
we saw in the Fourier coefficients of Figure [7|. These oscillating signs could be taken as an 
indication that an effective string representation of loops traversed multiple times requires 
some elementary two dimensional excitation of statistics different from that of bosons |m| . 
The simple effective string models one typically uses do not address loops traversed multiple 
numbers of times. 

Because of the importance of the large iV universality class to our program, we now 
turn to a more general description of what appears to be its fundamental representative 
family in terms of matrix models: 

Imagine n SU(N) matrices, independently and identically distributed according to a 
conjugation invariant measure dfi(Ui). The measure is peaked at XJi = 1 and is of the single 
trace type. As a function of each eigenvalue it decreases monotonically as one goes round 
the circle from 1 to —1, where it has a minimum. The continuum Wilson loop operators, 
W, one defines, would behave in the vicinity of the large ./V transition point at which their 
eigenvalue distribution opens a gap at —1, similarly to W = TJ™ U{. The parameters of the 
measure and n can be adjusted to various critical behaviors in the eigenvalue distribution 
of W at infinite N, and likely produce a variety of double-scaling limits corresponding to 
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Figure 9: Fit of the distributions to the Durhuus-Olesen distributions for four different sizes of 
Wilson loops, namely, lt c — 0.740, 0.660, 0.560, 0.503. The associated areas (k in the Durhuus-Olesen 
notation) that describe the continuous curves are given by k = 4.03, 2.30, 1.41, 1.15 respectively. 

large ./V fixed points of decreasing degrees of stability. The shape and scale of W map 
analytically into these parameters in the vicinity of the transition. Most plausibly, the 
most generic and stable transition is the relevant one. 

This random matrix model is exactly soluble, as we now sketch, leaving details for the 
future: 

One introduces n pairs of iV-component Grassmann variables, ipi,ipi and proves by 
simple recursion that 

/n 
H{d^dtPi)e z ^^-^^ u ^ +1 = det(z n + W) (4.7) 
i 

Above, we use the convention VVi+i = —ipi. 

Now, one can introduce auxiliary fields for the SU(N) invariant bilinears V'iV'i an d 
V'iV'i+i- Averaging over the Ui variables can be done independently and the remaining 
integral can be dealt with at infinite N by ordinary saddle point methods. 

Choosing an axial gauge on an infinite two dimensional lattice, it is easy to see that both 
on the lattice and in the continuum [18| the non-minimal Wilson loops of two dimensional 
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planar QCD are indeed described by the above universality class. Similarly, the Douglas- 
Kazakov critical point is also in this universality class, corresponding to the most stable 
critical point. 

Intuitively, the universality class can be visualized as a nonabelian generalization of 
the following abelian situation: Let there be given a smooth loop C and let S be a surface 
of the topology of a disk bounded by C. On the surface S define a two dimensional U(l) 
gauge field, with independently and identically distributed fluxes through small patches 
of surface that tile S. The phase of the parallel transport round C is then given by the 
sum of all these fluxes. The distribution of that sum will be generically governed by the 
central limit theorem. It is now obvious that we have a nonabelian generalization of this 
arrangement [ 31 1 , of the type studied by Voiculescu [32]. 



5. Summary and Discussion. 

We first give an overview of various ideas that have influenced us and to which we have 
not yet referred directly. Subsequently we proceed to a summary explaining what are the 
essential new points we have made in this paper. We finish with a short description of 
planned future work. 

5.1 Connections to other work. 

At the conceptual level, this paper has been influenced by the work of many others, in 
particular people in our own subfield of lattice field theory, even if this was not mentioned 
at any point in our presentation until now. While it is impossible to avoid omissions, we 
try here to draw attention to what stands out in our memory at this point. 

That an extra dimension might be needed for a local string dual of the field theoretical 
description clearly is an idea with many sources. The relationship between this extra 
dimension and a space of eigenvalues has been emphasized by many, for example in papers 
by A. Jevicki [33]. The essential point is that for a hermitian or unitary matrix the 



eigenvalues are restricted to a one dimensional line and repel each other; as such they 
make up dynamically, at infinite N, a continuous one dimensional degree of freedom. In 
the context of matrix models this produces an analogue of the extra dimension we have 
seen in higher dimensional setups. 

Using the lattice to make progress on the question of a possible string dual of non- 
abelian YM theory has preoccupied a sizable group of researchers. This new wave of lattice 
work on the large N limit has been started by M. Teper, and in various collaborations he 
has addressed many of the questions we address here [34]. In particular, the question of 



the cross-over has been studied in various contexts 1 35 ] . There is an interesting issue raised 
there suggesting a qualitative difference between three dimensional and four dimensional 
YM gauge theories with respect to observable dependent large N transitions. The imme- 
diate implication is that it would be interesting to repeat our work here in three Euclidean 
dimensions; it would be surprising to see three dimensions so different in this context from 
both two and four dimensions. 
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J. Kuti [gg] and collaborators have studied in great detail the cross-over from field 
theoretical to string-like behavior and have shown how various string states reorganize 
themselves as one goes from short to long QCD strings. The philosophy of matching various 
field theories with known or suspected string-like excitations onto an effective string theory 
at long distances underlies much of this work. This work has not yet been extended to 
large values of N, but this is recognized as a possible valid avenue for further steps in this 
program. 



We have been influenced significantly by papers by R. Brower [37]. On the one hand 
side he and collaborators have provided additional examples in which one can see the extra 
dimension being connected to standard scale decompositions of four dimensional processes 
of traditional interest. This interpretation is very natural when one thinks in terms of a 
warp factor in the context of the famous AdS/CFT correspondence in the conformal and 
related cases. It provides a most convincing resolution of the old puzzle of how a string 
theory can embody hard, field theoretic ultraviolet behavior. 

It has also been stressed by Brower that one should be more ambitious than merely 
connecting onto an effective long distance string theory, and, that the extra dimension 
could provide the means to construct a complete stringy dual to YM, valid at all distances, 
albeit too difficult for actual calculation at short distances. It is with this in mind that 
although our strings are fat, we maintain the zigzag symmetry, which would be a potential 
clue about the infinitely thin string dual. Also, with this in mind, we insisted on finding 
observables that have a continuum field theoretical semiclassical expansion at small sizes, 
where the effects of short distance fluctuations reveal themselves. 

The effective string description has been studied in various guises also by others p8|j . 
in addition to Polchinski and Strominger; the approaches focused both on open Dirichlet 
strings and on closed string wrapping round a large compact direction of space-time. 

Another related topic is the issue of "Casimir scaling" , in the context of the dependence 
of the string tension on the representation of the source and sink the SU(N) string connects. 
It is known that "Casimir scaling" is exact when applied to charges in the fundamental 
and in the adjoint representation at N = oo, so long as we take N — > oo first and let the 



minimal area spanned by the loop increase only after that 3£]. We wish to mention here 
a recent paper in the context of "Casimir scaling" at finite N where the Wilson loop 
eigenvalue distribution has been studied on the lattice for SU{2). 

Last but not least we wish to mention the old idea of preconfinement [[Tl|], in which 
one pictures a jet producing process in QCD as evolving by first creating (in the planar 
limit) finite energy clusters that already are color neutral and only subsequently hadronize. 
This preparation for hadronization is similar to the early signals of closing the gap in the 
eigenvalue spectrum, before the density becomes uniform up to corrections exponential 
in the minimal area spanned by the loop, where one has full confinement. As mentioned 
previously, one can have a large N transition without ever making it all the way to con- 
finement. 

5.2 Main new ideas in this work. 

Our main change in vantage was a focus on the universal nature of large N transition which 
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we hope to turn into a calculation that can take us through the cross-over which connects 
the regime in which field theory is the most convenient description to the regime in which 
string theory is assumed to eventually become a convenient description. 

We see the beginnings of such a calculational scheme and feel that were it successful, 
it would provide a concrete result on which to build in the search for a more elegant and 
fundamental dual representation of Yang Mills theory. 

We found ourselves naturally led into something that looks like an extra dimension 
and one that is intimately related to scaling at that. 

5.3 Plans for the future. 

The first problem for the future is a more thorough investigation of the critical regime of 
the large N transitions and the numerical identification of the right universality class. In 
parallel an effort will be made to make the semiclassical calculation based on instantons at 
the short distance end concrete. The double scaling limits in the larger universality class 
would need to be studied in detail; assuming that our guess for this universality class is 
right, it should be possible to construct those double scaling limits exactly. 

The connection onto a convenient effective string model at the strong coupling end is 
another direction of future activity; here, the main tool would be numerical at first, and we 
expect to draw from the experience of the many other workers on similar problems. Our 
main focus will be on keeping N large so the string coupling constant can be set to zero 
as is often done even for N = 3. 

Also, we feel it might be useful to try to face the question whether an effective string 
theory for planar QCD is indeed not different from an effective string theory in a model 
which has nothing to do with non-abelian gauge theory. For example, on the face of it one 
would expect the zigzag symmetry to play no role at leading order in the effective string 
theory, but maybe this is not a correct assumption. An analogy is the chiral Lagrangian 
describing pions in QCD, which indeed could come from a fundamental theory that has no 
gauge fields. But, the complete decoupling at all scales of all the mesons, well known to 
occur at infinite N, is something the effective Lagrangian, by itself, does not incorporate 
in a structural way; rather, many of its free parameters need to be set to zero. 
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